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Abstract 

Hydrogen arranges at dislocations in palladium to form nanoscale hydrides, changing the vibrational 
spectra. An ab initio hydrogen potential energy model versus Pd neighbor distances allows us to predict 
the vibrational excitations for H from absolute zero up to room temperature adjacent to a partial dislocation 
and with strain. Using the equilibrium distribution of hydrogen with temperature, we predict excitation 
spectra to explain new incoherent inelastic neutron-scattering measurements. At OK, dislocation cores trap 
H to form nanometer-sized hydrides, while increased temperature dissolves the hydrides and disperses H 
throughout bulk Pd. 

PACS numbers; 61.72.Yx, 61.72.Lk, 63.20.dk, 63.22.-m 
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I. INTRODUCTION 



The increasing needs for renewable energy — and issues of production, storage and transporta- 
tion of energy — motivates interest in hydrogen for energy storage.' At a fundamental level, open 
questions remain about how hydrogen acts in metals, despite a long legacy of study .■^ Palladium 
is an ideal metal to study hydrogen behavior due to the strong catalytic behavior of the Pd surface 
facilitating hydrogen adsorption, favorable T - pn^ thermodynamic properties, and that hydrogen 
acts as an ideal lattice gas in Pd.- Neutron-scattering characterization is useful,"^ in part to a scat- 
tering interaction mediated by neutron-nuclear properties and available incident neutron energies 
similar to those associated with lattice vibrations. Coherent inelastic neutron scattering gave the 
first phonon dispersion measurement of a metal hydride (Pd-H and -D).^ The hydrogen-dislocation 
trapping interaction in Pd has remained of significant interest over the last four decades because 
of the favorable Pd-H properties mentioned above and that Pd can be heavily deformed by hydride 
cycling across the miscibility gap.- 

Mobile solutes — substitutional and interstitial — arrange themselves in a crystal to minimize 
the free energy; with non-uniform strains, the arrangement reflects the energy changes from strain. 
For an edge dislocation, compressive and tensile strains produce areas that are depleted and en- 
hanced with solute concentration — a "Cottrell atmosphere."- Cottrell atmospheres produce time- 
dependent strengthening mechanisms like strain-aging in steels and the Portevin-Le Chatelier ef- 
fect in aluminum alloys,- and the rearrangement of hydrogen from dislocation strain fields aff'ects 
dislocation interactions.— The dislocation core — where the continuum description of the strain 
fields breaks down — provides the largest distortions in geometry and the attraction of solutes to 
this region is crucial for solute effects on strength.— — Tensile strain also lowers the vibrational 
excitation for H, and, in a dislocation core, broken symmetry splits the excitations.— The vibration 
of Pd next to H changes the local potential energy for each H atom, broadening the vibrational 
excitations. Additionally, the vibrational excitations of the light hydrogen atom are significantly 
changed by anharmonicity.-^^ We treat all of these effects: non-uniform hydrogen site occupancy 
due to strain and H-H interaction, quantum-thermal vibrational displacements for neighboring Pd, 
and the anharmonic potential energy to determine the causes of changes to the vibrational spec- 
tra with temperature. Experimentally, in situ inelastic neutron scattering averages over different 
H sites to give a direct measurement of H environment. We compare our ab initio treatment of 
hydrogen sites and anharmonic vibrational excitations with incoherent inelastic neutron-scattering 
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measurements to observe the formation and dissolution of nanoscale hydrides around dislocation 
cores in palladium. 



n. METHODS 

Incoherent inelastic neutron scattering (IINS) using the Filter Analyzer Neutron Spectrometer 
(FANS) at the NIST Center for Neutron Research measure the vibrational density of states of 
trapped hydrogen in polycrystalline Pd as a function of temperature. FANS scans the incident 
neutron energy and records the intensity that passes through a Be-Bi-graphite composite neutron 
filter. Sample preparation procedures and material are identical to Heuser etal.,^^ with ~ 100 grams 
of polycrystalline Pd sheet measured at 4K, lOOK, 200K, and 300K. Palladium sheet supplied by 
Alpha Aesar was cold-rolled in the as-received condition, and further deformed by cycling twice 
across the hydride miscibility gap.- It was held under vacuum at room temperature for several 
days and then annealed for 8 hours at ~400K to completely outgas the sample. The subsequent 
measured pressure reduction in a closed volume at room temperature using a portable hydrogen gas 
loading apparatus gives a total hydrogen concentration of 0.0013 [H]/[Pd], corresponding to a total 
hydrogen inventory of 1.3 mg. The IINS measurements were performed in an Al measurement 
can sealed with indium wire.— This can was isolated with an all metal vacuum valve, mounted 
to the FANS instrument, and cooled to 4K. Subsequent measurements were performed at lOOK, 
200K, and 300K. The sample was then outgassed at ~420K for ~48 hours completely remove 
all hydrogen. The zero-concentration background was measured from the out-gassed sample in 
the Al can at 4K, lOOK, 200K, and 300K. We also recorded fast neutron background with the 
sample in place and the detector bank blocked with Cd.^*^ The measured hydrogen vibrational 
density of states in Fig. [3] is the normalized net intensity after zero-concentration and fast neutron 
background subtractions. In addition, an energy-independent flat background attributed to multi- 
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phonon scattering was subtracted, as discussed in Ref. 

Density functional theory calculations for Pd-B^-"^ are performed with vasp^^ using a plane- 
wave basis with the projector augmented-wave (PAW) method^^ with potentials generated by 
Kresse.-^ The local-density approximation as parametrized by Perdew and Zunger— and a plane- 
wave kinetic-energy cutofi" of 250eV ensures accurate treatment of the potentials. The PAW po- 
tential for Pd treats the s- and rf-states as valence, and the H 5-state as valence. The restoring 
forces for H in Pd change by only 5% compared with a generalized gradient approximation, or 



including Pd 4p-states in the valence; our choice of the local-density approximation is computa- 
tionally efficient, and gives an a-Pd lattice constant of 3. 8528 A compared with the experimentally 
measured 3.87 18A. To compute the dynamical matrix for Pd, and to relax H at the octahedral site 
in Qf-Pd, we use a 4 x 4 x 4 simple-cubic supercell of 256 atoms, with a 6 x 6 x 6 k-point mesh; 
while the dislocation geometry with 382 atoms uses a 1 x 1 x 8 k-point mesh. For the PdHo.ga 
hydride force-constant calculation, a 3 x 3 x 3 simple cubic cell (108 Pd atoms, 68 H atoms) with 

D 

displacements of 0.01 A for H and Pd atoms and a 8 x 8 x 8 k-point mesh. The electron states are 
occupied using a Methfessel-Paxton smearing of 0.25eV. For the H octahedral site in a-Pd and the 
partial dislocation core, atom positions are relaxed using conjugate gradient until the forces are 
less than 5meV/A. 

Dislocations produce a distribution of interstitial site strains; to compute the density of strain 
sites available for hydrogen, we consider a simplified model for the distribution of dislocations 
throughout the crystal. We take the dislocation density p^isi as given by cylinders of radius R = 



1 / V^rpdisi with an edge dislocation at the center; we assume that the strain in each cylinder is due 

only to the single edge dislocation at the center. The volumetric strain r away from the dislocation 

core and with angle 6 to the slip plane is 

^ l-2v . ^ . ^ b 

£ = -- — • -; smO = - smO ■ - — (1) 

2nr I - V 4nr 

for a Poisson's ratio v = 1/3, and where b = 0.298nm is the Pd Burgers vector. This equation 

becomes invalid for small r; we truncate the expression in the "core" of the dislocation. We can 

estimate the size of the core by considering the maximum strain of +5% at the partial core from 

Ref. bj; then. 

The line vector of an edge dislocation is t = V6ao/2 with Burgers vector b = ao/ V2, and so the 
core has a volume of r^t = 3aQ/2 = 6(ao/4); hence, there are 6 sites per dislocation line inside 
this radius. We assign half the maximum strain of +5% and half the minimum strain of -5% 
corresponding to opposite sides of the partial cores. Previous ab initio calculations of the core 
give a trapping energy of 0.164eV with a 5% strain; the trapping energy matches the decrease 
in hydrogen energy from a 5% increase in volume — we then model the binding energy for H as 
linear in the site strain s: -0.164eV(e/0.05). 

With these definitions, we compute the density of strain sites by integrating over our cylin- 
der cross-section from out to R. We consider the 6 core sites (3 attractive and 3 repulsive) 



separate from this continuum calculation. 
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where the delta-function integral is calculated by rewriting the delta function in terms of the two 
roots 6q = sirT^sAnr/b). To simplify the expression, we define two strains: the maximum site 
strain si = bliAnr^), and the maximum strain at the cylinder edge eq = b/(4nR). Then, 
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These two expressions can be written in terms of the ratio t] = eq/ei < I as 



n(E) = 
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The general scaling n ~ |e| ^, similar to Kirchheim.^'* If we integrate this density of states over all 
strains, we have 



which accounts for the "missing" core states, which are a fraction = r^/R^ of all possible sites. 
We add back the core sites that make up 6r^pdisi of all possible sites; half have tensile strain +si, 
and the other half have compressive strain -si. In our sample, the dislocation density is pdisi = 
lO'^cm"^, so R = 1/ y;rpdisi = = 19nm, the maximum site strain is Si = bl{Anrc) = 0.05, 

and the maximum strain at the cylinder edge is eq = bl{AnR) = 1.25 x 10^^, with a ratio of 
T] = so/si = 0.025, and with a core occupancy of 6r^pdisi = 2 • 0.576 x 10"^. 

The thermodynamics of hydrogen in Pd requires considering not just the site strain from a dis- 
location, but also from neighboring hydrogen atoms. The site adjacent to a hydrogen interstitial 
in Pd experiences strain due to the occupancy of the hydrogen site; this strain, in term, affects 
the site energy. In a 256-atom Pd supercell calculation of a hydrogen interstitial, the relaxation 
neighboring the hydrogen interstitial site is expanded by Ae = 6.864 x 10"^; this produces a low- 
ered site energy of approximately IS.E = -23meV. It should be noted that this is purely classical 
approximation — it ignores not only electronic structure effects, but zero-point displacement of the 
two hydrogen atoms. However, it should give the correct order of magnitude for the strength of 
interaction, and it suggests a propensity for ordering on the hydrogen sublattice. 

To account for the weak H-H binding on the hydrogen distribution and site occupancy, we 
consider a simple self-consistent mean-field model. A site with energy E (or, alternately, strain e) 
will be shifted by AE if any of its neighbors are occupied, and unshifted if all are unoccupied. We 
will ignore spatial variations in the local site occupancy, and so approximate the probability of each 
neighboring site being occupied with the site occupancy /. As there are twelve possible nearest- 
neighbor sites in the FCC hydrogen sublattice, the fraction of sites where all twelve neighbors are 
unoccupied is (1 -/)^ with = 12; hence, each site now has two possible energy levels: a fraction 
(1 ~ f)^ with energy E and a fraction 1 - (1 - f)^ with energy E + AE. To be in equilibrium, 
these sites have occupancies of /o = (exp(J3(/j. - E)) -\- 1)'^ and /i = (exp(J3(ju - E - AE)) + 1)"\ 
respectively. Thus, the occupancy of a site satisfies the self-consistent equation 



This equation is solved for / at each site given its energy E, and the chemical potential ju; the 
occupancy is integrated over the density of sites to determine the total concentration of hydrogen. 
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(8) 



/=/i+(i-/y(/o-/i) 



(9) 



Eqn.|9]can be solved approximately (to 10 by making a quadratic approximation around f ~ f\ 
to / = g{f). Defining the function and its first two derivatives at /i, 

^^°' = /i + (i-/iA/o-/i) 

g^') = -N{l-f,f-\fo-f,) (10) 
g^^y=N{N-m-fyf-\fo-h) 
the quadratic approximate self-consistent solution is 



/ = 2 



1/2 



(11) 



This self-consistent mean-field model accounts for the hydrogen-hydrogen attraction, and the pri- 
mary eff"ect is at low (but above zero) temperature where the ordering competes with entropy; it 
produces somewhat higher hydrogen occupancies than would be expected without any H-H in- 
teraction. This approximate thermodynamic model is not accurate when the hydrogen occupancy 
becomes large; for example, it does not account for the formation of PdHo.es before the formation 
of PdH. 

Fig. [H shows the formation of Cottrell atmosphere at low temperatures and dissolution near 
room temperature, including the difference between integrated occupancies assuming ls.E = 
and AE = 23meV. Qualitatively, assuming AE = shows similar behavior to AE = 23meV, 
with dissolution of the nanoscale hydride between 200K and 300K. The primary eff"ect of the H-H 
binding is to maintain a slightly higher hydrogen concentration in the dislocation cores. Fig. [T] 
shows the integration of site occupancy, starting from the core; the derivative with strain gives the 
fraction of H at a specific strain. At OK and lOOK the core is fully occupied; hence, the integrated 
occupancy starts at 0.576 x 10 At 200K the core is 96% occupied, falling to 54% occupancy 
at BOOK. As temperature rises, lower strain sites have an increased occupancy due to entropy, and 
sites near the core are less populated — the "dissolution" of the Cottrell atmosphere, though the 
core still has hydrogen. The fractional occupancy of sites near zero strain decays exponentially, 
but as the number of sites is growing as |e|~^ most of the hydrogen is well dispersed at higher 
temperatures. 

Prediction of vibrational excitations for hydrogen requires sampling of diff"erent Pd displace- 
ments neighboring the H atom to determine the potential energy. Hydrogen is surrounded by 6 
Pd neighbors at |(100). These six neighbors are displaced according to the thermal occupation 
of phonons, including the quantum-mechanical zero-point motion. The displacements provide an 
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radius = h/4ne 
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FIG. 1. (Color online) Integrated occupancy of hydrogen around dislocations in Pd with temperature for 
AE = 23meV (solid) and AE = (dashed). The integrated occupied density of sites goes from the most 
favored sites (dislocation cores) through the range of volumetric strain around the dislocation core; all 
hydrogen solutes are accounted for at the saturation concentration of .xh = 1.3 x 10~^. The occupancy 
follows a Fermi function for AE = 0, and the effect of H-H coupling is to maintain the nanoscale hydride 
to slightly higher temperatures. At OK, the Cottrell atmosphere has a sharp boundary at r = 4^ = 1.9 A. At 
lOOK, the atmosphere shows only small spreading away from the core, while at 200K there is an increasing 
occupancy for H at strain. At 300K, the atmosphere is dissolving, with decreased occupancy in the core 
as well as around the dislocation. 



important broadening of the hydrogen vibrational excitation spectra, as the light hydrogen atom 
evolves in a Bom-Oppenheimer-like manner (valid as Mh « 10"^Mpd), sampling the local po- 
tential energy from the neighboring Pd. To compute a density of excitation energies for the H 
atom, we need to sample the possible displacements for neighboring atoms at a temperature T. 
For the highest frequency excitation of Pd, 8THz {Hco = 33meV), Xq/ V2 = 0.025A; at 300K, 
X = 0.033A. The Gaussian distribution of displacements for a harmonic oscillator (see Appendix) 
provides the basis for random sampling displacements for Pd atoms from independent Gaussians 
of width x{co„{q), T) for each phonon mode a»„(^ in the Brillouin zone. Let D{R) be the 3x3 
force-constant matrix between an atom at and ^; moreover, let U{^) be the displacement vector 
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for an atom at R. Then, the Fourier transforms of D and u are 

R 

for a bulk system of A'^ atoms. The inverse Fourier transforms are 

D{R) = ^Yj^icOe-'^l-^ 



(12) 



(13) 



where we have used the fact that there are also q-points in the Brillouin zone summation. Note 
also that, 

^ = J] . (14) 

it 4 

Then, the displacements can be written as the sum of three Gaussian distributed random 
variables Uniq), multiplied by the corresponding width x((x>n(q), T) and normalized eigenvector of 
D{q), Uniq)- In reciprocal space, the sampled displacement is 



n 



coth(na)n{q)/2kBT) 



1/2 



(15) 



ImtOniq) 

The final step is to inverse Fourier transform all of the displacements, and to remove the center- 
of-mass shift for the the six neighbors surrounding the H atom at {r\. In the sum over the discrete 
q in the Brillouin zone, the weight of each point w{q) = 1 /N, so 



Au(r) = ^ a„(q)Un(q) 



nq 



w{q)fi 



coth {fia)n{q)/2kBT) 



1/2 



.|cos(^.^_i^cos(<^.^)| (16) 



This requires 3A^ - 3 random Gaussian variables a„(q) to produce one sample of displacements for 
Pd atoms neighboring the hydrogen atom at a temperature T. 

The force-constants for Pd come from ab initio via a direct-force technique- with a 4 x 4 x 4 
simple-cubic supercell; this reproduces the elastic constants and phonons within 5%. We use a 
discrete 16x16x16 Monkhort-Pack mesh^ of ^-points the Brillouin zone. With 40,000 dis- 
placements for each temperature (OK to 300K), in the dislocation core and strains from -1-0.05 to 
-0.01 in 0.01 increments, we compute vibrational excitations for H in Pd. Given the H potential 
energy, we solve the Schrodinger equation numerically. For each Pd displaced environment, we 



find the minimum energy position for H, and expand the potential as a fourth-order polynomial 
in H displacement, and compute the three lowest-lying excitations using a Hermite-polynomial 
basis.— This gives 120,000 excitation energies, binned into ImeV bins. Thus, we predict vibra- 
tional density of states for H in a dislocation core, and at strains from -1-0.05 to -0.01 at OK, lOOK, 
200K, and 300K. 

To efficiently describe the energy landscape for a hydrogen atom in a variety of interstitial 
sites — including small displacements of Pd due to quantum-thermal vibrations — we optimize an 
embedded-atom method-like potential for H based on its distance to six neighboring sites. The 
embedded-atom method^^— can work well for describing the energy of atoms in metallic systems: 
neighboring atoms have overlapping charge densities at a site, and atoms experience an "embed- 
ding energy" due to that local environment. As we are interested in describing H accurately for 
a small range of environments, we define a potential based on similar ideas, but make the fitting 
parameters as linear as possible so that overfitting can be easily identified, and good transferability 
achieved. From previous calculations,'"^ we have a large amount of force-displacement data for H 
in diff"erent environments (58 displacements in the dislocation core, 40 displacements in unstrained 
Pd, and 32 displacements in -1-5% strained Pd). This fitting database gives sufficient coverage that 
our potential will be used to interpolate rather than extrapolate. The general form of the total 
energy in terms of the H-Pd distances is 

D 

d=2 in 



where 




m 



where D and C determine the polynomial order of the embedding energy U (p) and the pair po- 
tential (pir); besides the coefficients Ud and cpc, there is the parameter a which determines decay 
length of the density. This means that the energies (and forces) are linear in all parameters except 
a; we can easily optimize the parameters by solving for Ud and 0^ for a given a with the smallest 
mean-squared error in the forces (weighted by the force magnitude). Hence, for any choice of D 
and C, we can find optimal parameters to accurately reproduce the DFT forces. To optimize the 
choice of D and C, we computed the leave-one-out cross-validation score (CVS) for each optimal 
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FIG. 2. (Color online) Calculated vibrational density of states for hydrogen in Pd with temperature. Increas- 
ing temperature produces larger displacements of Pd beyond the zero-point motion at OK; this increases the 
spread in the vibrational excitations. The central peaks for the three sites are temperature independent. Peak 
broadening smears the low and high excitations in the dislocation core at room temperature. 

set of parameters; D = 2 and C = 5 had the lowest CVS. This fit {Ey^ in eV, in A), 

(18) 

+ Yj[~ 131.94r„, + 119.41r,^, - 54.073r,^, + 12.1883r^, - 1.09167ry 

m 

had no error larger than 10% in any of the forces, and reproduced the H excitation spectra of the 
direct DFT calculation to within 2meV. As r < 2A, the contribution of the higher order polynomial 
coefficients is decreasing to larger orders. 



III. RESULTS 

Fig. [2] shows the predicted vibrational density of states for hydrogen at equilibrium zero strain, 
a 5% expanded site, and in the partial core. Increasing temperature broadens the excitation spectra 
with increased vibration of neighboring Pd atoms. There is no shift in the peak position with 
temperature due to Pd vibration, but only from strains. The dislocation core environment breaks 
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cubic symmetry, giving three peaks below 120meV.— Temperature widens the peaks above and 
below 78meV on each side of the central peak at room temperature. Hence, despite dislocation 
core occupancy at room temperature, it is difficult to experimentally identify H in the dislocation 
core except at low temperatures. 

Fig. [3] shows the predicted vibrational spectra for 0.13at.% H in Pd as a function of temper- 
ature, and the comparison with inelastic neutron scattering measurements. Combining the site- 
occupancy data from Fig. [T] with the predicted vibrational spectra in Fig. |2l we predict the expected 
measured vibrational spectra with temperature. To compare with the experimental measurements, 
we scale all of our peak heights to be equal, scale intensity by 1 / ^Jhv to produce a scattering 
cross-section under the condition of variable incident energy and fixed final energy (as is the case 
for the measured TINS spectra reported here), and scale energy by 7/8. The latter scaling corre- 
sponds to a needed softening of the DFT calculations of vibrational spectra for H in Pd compared 
with experimental measurements; the overestimation of vibrational excitation is independent of 
exchange-correlation potential and treatment of H and Pd ionic cores^ and is consistent with ear- 
lier fuUy-anharmonic calculations of isolated hydrogen in Pd.— The experimentally measured line 
shape is in good agreement with the prediction of scattering at room temperature, but the shapes 
begin to deviate as temperature is lowered. 

Lowering temperature forms a Cottrell atmosphere and the predicted scattering cross-section 
shifts and narrows; the shift in peak energy agrees with the experimental measurements, but the 
peak narrowing does not. At 300K, hydrogen is primarily in low strain environments, and has a 
peak widened primarily by vibration of Pd neighbors. As temperature is lowered, ab initio calcu- 
lations predicts a shift of the peak to lower frequencies as higher strain sites and the dislocation 
core is preferentially occupied; this matches the experimental measurement as well. However, the 
ab initio calculations predict a narrowing of spectra; this narrowing is due to the smaller displace- 
ments of Pd neighbors producing less random distortion of the potential energy. As the Cottrell 
atmosphere forms, the local hydrogen concentration near the dislocation core is very high, form- 
ing hydride phases in nanoscale cylinders. This corresponds well with recent small-angle neutron 
scattering measurements at low temperatures and hydrogen concentrations in deformed single- 
crystal Pd.— The vibrational spectrum of jS-PdH is wider due to H-H interactions;— this dispersion 
is lacking in the ab initio calculations due to the difficulty of predicting fully anharmonic dis- 
persion relations. We have computed the harmonic bulk PdHo vibrational density-of-states that 
includes dispersion, but lacks anharmonicity; the comparison with the IINS signal from OK to 
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FIG. 3. (Color online) The predicted vibrational density of states and inelastic neutron scattering intensity 
for 0.13at.% H in Pd as a function of temperature. The temperature determines both the occupancy of states 
for H (c.f. Fig.[Tll and the vibrational spectra for all states (c.f. Fig.O; taken together, we predict the density 
of states in the top figure. To compare with IINS measurements, we scale intensity by 1/ ^^hv, equalize 
amplitudes, and scale energy by 7/8 (DFT/experimental discrepancy). The agreement in line shape at 300K 
confirms that the main cause of peak broadening is Pd vibration. At lower temperatures, the formation of a 
Cottrell atmosphere creates nanoscale regions with high hydrogen concentration. The scattering signal from 
yS-PdH has a width similar to the experimentally measured spectrum at OK;— the difference from the ab initio 
prediction is due to the dispersion of a hydride which is missing in our calculation of isolated hydrogen 
vibrations. The computed PdHo.63 spectra (dashed lines) has dispersion but is a harmonic approximation 
for hydrogen. The signal change can estimate the fraction of nanoscale hydrides at dislocation cores. 
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200K strongly backs up the presence of hydride. Fitting the experimental intensity to a linear 
combination of the two predicted intensities suggests all hydrogen is in hydride and none is free 
at OK; a 9:1 ratio at lOOK; a 9:4 ratio at 200K; and dissolution of the hydride at 300K. Hence, 

we conclude that the Cottrell atmosphere is forming of nanoscale hydride particles near disloca- 
tion cores, despite the low total hydrogen concentration in the sample, to explain the changes in 
vibrational spectra. 

IV. CONCLUSION 

Combining the experimental measurement of hydrogen vibrational spectra with ab initio cal- 
culations of vibrational spectra with temperature, we can identify the formation of Cottrell atmo- 
sphere leading to nanoscale hydride precipitates at dislocation cores. By separating the sources 
of spectral broadening — dispersion in hydrides at low temperatures, and thermal broadening from 
Pd vibration of neighbors — and the causes of a peak shift, we have in situ characterization of the 
hydrogen environment evolution with temperature. 
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Appendix A: Harmonic displacement distribution at finite temperature 

For an isolated harmonic oscillator of mass M and natural frequency oj, we want to determine 
the probability distribution of displacements x from equilibrium. The state energies are En = 
fia){n + 1 12), and so the probability of being in state n at temperature T (fi = (ksTy^) is 




(Al) 



The wavefunctions are 



ilr„(x) = (2"n\r"\nxlr"'e-'^'^^oH„(x/xo) 



(A2) 
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for natural length xq = ylWjmco, and Hermite polynomial Then the probability distribution of 
displacement x is 



n=0 



e-^'I'^lHl^xIxo) 

^^,) 2"n\ ^/nxo 



1 _ 1 

1 



exp - 



1 - g-^'*'^ 
1 + e-f^''"-' xl 



(A3) 



exp 



V2^Jc(a;,r) \ 2x\io,T) 



where 



x((t>, T) 



n 



2ma> 



coth 



1/2 



(A4) 



is the thermal Gaussian width; the simplification is possible by using Mehler's Hermite polynomial 
formular^ 



Hn{x)H„{y) (wY 



w^r^'^ exp 



2xjw - (^x^ + y^) 

1 -W2 



In the low temperature limit, x ~ xq/ y/2 as expected from zero-point motion; and in the high 
temperature limit, x « (kBT /mo/y^^, as expected from the equipartition theorem. 
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